Relativistic fluctuating hydrodynamics with memory functions and colored noises 
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Relativistic dissipative hydrodynamics including hydrodynamic fluctuations is formulated by 
putting an emphasis on non-linearity and causality. As a consequence of causality, dissipative 
currents become dynamical variables and noises appeared in an integral form of constitutive equa- 
tions should be colored ones from fluctuation-dissipation relations. Nevertheless noises turn out to be 
white ones in its differential form when noises are assumed to be Gaussian. The obtained differential 
equations are very useful in numerical implementation of relativistic fluctuating hydrodynamics. 
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Introduction — Relativistic hydrodynamics has been 
widely used so far in various fields such as cosmology, 
astro- and nuclear physics to describe space-time evolu- 
tion of thermodynamic variables and flow phenomena. In 
the physics of relativistic heavy ion collisions, relativistic 
hydrodynamic models play a crucial role in drawing a re- 
markable conclusion that the quark gluon plasma (QGP) 
[1], which is novel deconfined nuclear matter and filled 
the early universe, is created very quickly after the col- 
lisions [2] , records the highest man-made temperature of 
the order of lO^^K [3, 4], and, most interestingly, behaves 
like a perfect fluid [5] . A new paradigm of strongly cou- 
pled QGP was established through a vast body of efforts 
of both the experimental [6-9] and theoretical analyses 
[10-13]. The next task in this field is to precisely extract 
transport coefficients of the QGP such as shear and bulk 
viscosities through comparison of results from relativistic 
hydrodynamic simulations of the QGP with experimental 
data intimately related with collective flow phenomena. 

In ordinary hydrodynamic framework, thermodynam- 
ics is assumed to be applicable in a local rest frame at 
each space-time point. This implicitly means that there 
exists an intermediate scale at which the size of a fluid el- 
ement can be characterized. At this scale, a microscopic 
characteristic scale should be small enough to apply ther- 
modynamics in a fluid element, while macroscopic char- 
acteristic scale should be much larger than that scale. In 
order to obtain a set of equations of motion and consti- 
tutive equations, one can introduce the Knudsen num- 
ber, K = l/L, where I is the microscopic scale and L is 
the macroscopic scale, and expand macroscopic quanti- 
ties in a series of the small parameter K . In this way, hy- 
drodynamics has been applied to macroscopic system to 
describe slow and long-wave-length-limit dynamics. Hy- 
drodynamics is applicable when separation of these three 
scales can be justified. 

Since thermodynamics can be derived by taking an 
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ensemble average of many canonical systems in statis- 
tical physics, this concept lies implicitly also in hydro- 
dynamics. In the physics of relativistic heavy ion col- 
lisions, hydrodynamic simulations had been performed 
from smooth initial fields [14, 15]. The resultant dynam- 
ics is regarded as an average behavior of many similar 
events with almost the same impact parameter. Recently, 
event-by-event hydrodynamic simulations have been per- 
formed extensively to describe experimentally measured 
momentum anisotropy of hadrons [16]. When the QGP 
is created ~ 1 fm/c just after the collisions, causality 
does not allow the coarse-grained scale to be larger than 
^ 1 fm. Consequently the system responds hydrodynam- 
ically to the granular structure originating from fluctuat- 
ing configuration of nucleons inside colliding nuclei. How- 
ever, it is far from trivial whether thermodynamic and 
hydrodynamic concepts within its current form can be 
compatible with event-by-event description of relativis- 
tic heavy ion collisions. Since information of the system 
is reduced in the coarse-grained process, thermodynamic 
variables fluctuate from event to event around average 
values, which is nothing but thermal fluctuation. There- 
fore the thermal fluctuation should be taken into account 
towards a consistent description of dynamics of thermo- 
dynamic fields on an event-by-event basis. 

The thermal fluctuations which arise in hydrodynamic 
evolution independently for each spatial and temporal 
points are called hydrodynamic fluctuations. A non- 
relativistic theory of hydrodynamic fluctuation was initi- 
ated by Landau and Lifshitz long time ago [17, 18]. They 
extended the Navier-Stokes equations by implementing 
stochastic fluxes into constitutive equations based on 
fluctuation-dissipation relations. The theory treats linear 
fluctuations around equilibrium state. There have been 
several attempts to study non-linear hydrodynamic fluc- 
tuation within the framework of Fokker-Planck equation 
[19] or that of Langevin equation [20]. Non-linearities are 
of particular importance for fluids undergoing, e.g., phase 
transition, nucleation, and instabilities since these can 
amplify the effects of fluctuation exponentially [21, 22]. 

The effect of thermal fluctuation on transport co- 
efficients in the relativistic system was discussed in 



Refs. [23, 24]. An extension of linear hydrodynamic fluc- 
tuation in the Landau-Lifshitz theory [17, 18] to the rel- 
ativistic hydrodynamics was discussed in Refs. [25-27]. 
Hydrodynamic fluctuations on top of one-dimensionally 
expanding background field together with its implication 
to the physics of relativistic heavy ion collisions were also 
investigated [26, 27]. 

On the other hand, a naive relativistic extension of the 
Navier-Stokes equation [17, 28] has a problem on causal- 
ity [29] . The second order correction terms to the entropy 
current play an essential role for the theory to obey the 
causality [30, 31]. Thus the non-linear version of rela- 
tivistic fluctuating hydrodynamics being consistent with 
causality is demanded for the purpose of event-by-event 
hydrodynamic simulations of, e.g., relativistic heavy ion 
collisions. In this Letter, we formulate relativistic fluctu- 
ating hydrodynamics by putting emphasis on the impor- 
tance of causality, its consequence in the property of hy- 
drodynamic fluctuation, and non-linearity of the result- 
ing equations of motion. In the following, the Minkowski 
metric is g^^ = diag(+, — , — , — ) and the natural unit 
/i = c = /cs = 1 is employed. 

Memory function and fluctuation- dissipation rela- 
tion — Hydrodynamic fluctuation is deviation of the 
dissipative currents, such as bulk pressure and shear 
stresses, from their average values determined by con- 
ventional constitutive equations. Hence the effect of hy- 
drodynamic fluctuation appears only quadratic (and in 
general higher order) terms of hydrodynamic fields in 
the linearized fluctuating hydrodynamics. However it 
is far from trivial what happens when hydrodynamic 
non-linear equations are solved by taking account of 
hydrodynamic fluctuations. The hydrodynamic fluctu- 
ations are treated as stochastic processes in the hydro- 
dynamic equations. Thus the equations are no longer de- 
terministic ones, but stochastic differential equations like 
Langevin equations. Here we consider the Gaussian noise 
as the hydrodynamic fluctuations. The power spectrum 
of the hydrodynamic fluctuations is related to transport 
properties through the fluctuation-dissipation relation. 

To determine the power spectrum of the hydrody- 
namic fluctuations in relativistic dissipative hydrody- 
namics with flnite relaxation times, we can write down 
the constitutive equations explicitly with respect to the 
dissipative currents 



n 



nu 



d'^x' Gn{x - x')e{x') + 5Ii, 



(1) 



x°>x'° 



d^x'GAx-x'r^f'{d(^up^l,) + 5^^^\ 



[ d^x'G^J{x-x'f{T\/c^\^')+SL 

JxO>x'0 



(2) 



(3) 



where u^ is the four- velocity in the Landau frame [17], 
and T and ^i are temperature and chemical potential 
of the i-th conserved current (z = 1, . . . , n), respectively. 



The dissipative currents H, tt'^", and v^ are bulk pres- 
sure, shear stress tensor, and diffusion of the i-th con- 
served charge, respectively. The corresponding thermo- 
dynamic forces are 9 = d\u^, d/^uas = AaB-wsd'^u^ , and 
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TAc.fjd'^^, where A"'^ = g""^ - u^'u^ and 

jectors to each tensor component of the dissipative cur- 
rents. As above, the dissipative currents can be gener- 
ally written in the form of the average behavior of the 
response from the thermodynamic force plus its hydro- 
dynamic fluctuation SH, (^tt^", and Si'l^ . The average 
behavior can be expressed as a convolution of the ther- 
modynamic force and the retarded Green's function. The 
Green's function contains the information about trans- 
port properties of the system such as the shear viscosity 
rj, the bulk viscosity C, and the charge conductivities Hij. 
The function also contains the memory effect of the sys- 
tem such as relaxation time tr. Thus it is also called 
memory function. 

Using the fluctuation-dissipation theorem, we can ob- 
tain the power spectrum or the two point correlation of 
the hydrodynamic fluctuations as 

{Sn(x)SU{x')) = TG*u{x - x'), (4) 

{5n^"'ix)5T:"^{x')) = TGl{x - x'Y'"''^ (5) 

{5y^{x)5y-{x'))=TG:,{x-xr''. (6) 

Here the memory functions are extended for the do- 
main of x'^ < x^ to be even functions, e.g., G*Jx)'^°' = 
G*,(-a;)"^ = G.,(x)^" + G,,(-x)"^ 

In the case of the first-order dissipative hydro- 
dynamics with hydrodynamic fiuctuation, G^{x) — 
2C<5(4)(a;), Glix)''''"^ = 4:r]S'-^\x)A''''°'^ , and G*^(a;)^" = 
—2KijS^^'{x)A^°'. These memory functions arc local in 
time and have no memory effects: The correlation of the 
fluctuation vanishes in an infinitesimally small time du- 
ration. Thus the fluctuations exhibit white noises in the 
first-order dissipative hydrodynamics. 

Note that, in a more general case which respects the 
causality, the memory function should vanish with two 
spatially separate points, i.e., G{x — x') — where (x — 
x'Y < 0, while the delta functions in the first-order case 
are symmetric in any direction oi x — x' . 

A second-order dissipative hydrodynamics with hy- 
drodynamic fluctuations — As mentioned before, the 
first-order relativistic dissipative hydrodynamics violates 
causality [29]. To respect the causality in the relativistic 
theories, second or higher order dissipative hydrodynam- 
ics with relaxation effects should be considered. Here we 
consider a simple case of the second-order constitutive 
equations [30, 31] 
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where D = u^da is the time derivative in the Landau 
frame. The memory functions can be obtained by solving 



the constitutive equations for dissipative currents [32] 
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where a^ — (r(a;),cr(x)) are the proper time and the co- 
nioving coordinates in the Landau frame, respectively, 
and Texp(— J , dT"T~f, \t") is the time-ordered exponen- 
tial. The function 6* W(a;-a;') = |||;| j(3)(cr - cr')e(r - 
t') is defined as the common part of the memory func- 
tions. The tensors A(rf;Ti)^''"'' and A(Tf;Ti)^" are the 
time-by-time projection into the tensor components of 
the dissipative currents and defined as follows: 

A(Tf;ri)^''„^= hm A(rf)^"^oo^o 
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This non-trivial tensor structure comes from the projec- 
tions A^'^Q/3 and A^^ in the first term in the left hand 
sides of constitutive equations (8) and (9). This time- 
by-time projection is required to make tt^'' and i^l* con- 
strained to each tensor space during time evolution. 

In these memory functions (10)-(12), the time correla- 
tion has a form of exponential relaxation. In general, the 
memory function of causal relativistic dissipative hydro- 
dynamics has non-vanishing values for different times. 
Consequently, the time correlation of the fluctuations 
TG{x) does not vanish between different times, i.e., the 
fluctuations exhibit colored noises. For the bulk pressure, 
colored noises are manifested in the Fourier space as 
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Here (w, k) are the conjugate variables for (t, cr), respec- 
tively. For shear stress tensor and diffusion, the explicit 
formulas are much more involved due to time-by-time 
projection. Nevertheless, the corresponding power spec- 
tra in the Fourier space depend on uj similar to the bulk 
pressure case. 

Implementation in numerical simulations — There are 
two different strategies to implement hydrodynamic fluc- 
tuations in numerical simulations. One is the way to 



memorize the past fluid fields and to directly solve the 
constitutive equations with integral. The other is to go 
back to the differential form of constitutive equation such 
as Eqs. (7)-(9) with the hydrodynamic fluctuations. 

Since constitutive equations containing memory func- 
tions are non-Markovian, solving the equations requires 
the information about the all past fluid fields. So, in nu- 
merical simulations, one has to memorize the past fields 
of the all steps. Nevertheless one can introduce some cut- 
off step number of the past field to be memorized since 
the memory function relaxes to zero in finite time and 
the system should "forget" the information about past 
near the equilibrium. 

In the case of the fluctuating hydrodynamics with 
second-order dissipative terms discussed above, the con- 
stitutive equations become 
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It is difficult to evaluate exactly the time-ordered expo- 
nential and time-by-time projections, A(Tf; ri)^'^"'' and 
A(Tf ; Ti)''". So one can use the approximated expressions 
for TT^"^ and i^f 
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These equations are consistent with the constitutive 
equations of the differential form, (7)-(9), up to the sec- 
ond order in thermodynamic forces. 

A more practical way of performing numerical calcu- 
lations is to employ differential forms of the constitutive 
equations. For the simplest case of the second-order dissi- 
pative hydrodynamics, the constitutive equations includ- 
ing hydrodynamic fluctuations become 






(21) 
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Those equations are solved by regarding the dissipative 
currents as the dynamical variables in addition to the 
other fluid fields. This is the consequence of causality and 
the essential difference between the first order and the 



second order theories. It is shown that the noise terms in 
the constitutive equations in the differential form above 
become white ones even though the hydrodynamic fluc- 
tuation itself is colored one: 

(Cn(x)en(x')) = 2TC<5W(x-x'), (24) 

{^r{^)C^{^')) = 4rry<5(4)(x - a;')A^"^"^ (25) 

(er (x)e;(x')) = 2Tk,,5^^\x - x')A^". (26) 

Further more, it is also shown that, if the hydrodynamic 
fluctuations are Gaussian noises, the differential form of 
constitutive equations should have the following general 
form in the Fourier space: 
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Here H' is a dissipative current and nF is the first- 
order term of the thermodynamic force in the constitutive 
equation. Also, the noise term ^^^ /, (= [icjAfe + l](5n^ j^) 
becomes Gaussian white noise: 
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Therefore it is much more convenient to solve the equa- 
tions in the differential form in numerical calculations 
rather than the ones in its integral form if and only if 
a general assumption that noises are Gaussian can be 
made. The proof will be discussed elsewhere in the fu- 
ture publication [33]. 

Conclusions — To discuss about fluctuating hydrody- 
namics, which is the hydrodynamics with dissipation and 
the corresponding hydrodynamic fluctuations, in rela- 
tivistic system, one should consider the memory effects 
such as relaxation times to respect the causality. The 
constitutive equations can be written down in the form of 



the response of thermodynamic forces in the present and 
past, and the response functions are called memory func- 
tions. This is the generalization of the conventional first- 
order or second-order constitutive equations and they are 
the special case with certain memory functions. The 
power spectrum of the hydrodynamic fluctuations is re- 
lated to the memory function by fluctuation-dissipation 
relation. It should be noted that the memory function 
can be obtained in principle from the first principle calcu- 
lations such as lattice quantum chromodynamics [34, 35] 
through correlation function of energy momentum ten- 
sor in equilibrium. In the case of the first-order dissi- 
pative hydrodynamics, the fiuctuations are white noises 
while in the case of the second-order dissipative hydro- 
dynamics respecting causality, the fluctuations should be 
colored. If the hydrodynamic fluctuations are Gaussian 
noises, the differential form of the constitutive equations 
within the second-order relativistic dissipative hydrody- 
namics should have the specific structure, and the corre- 
sponding noises turn out to be white ones. 

Instead of conventional relativistic dissipative hydro- 
dynamic equations, the equations obtained in this Let- 
ter can be employed to simulate more realistically, e.g., 
space-time evolution of the QGP in relativistic heavy ion 
collisions on an event-by-event basis. Numerical imple- 
mentation of these equations together with the effects 
of hydrodynamic fluctuations on observables will be dis- 
cussed elsewhere [33]. 
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